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ON  THE  SECOND-ORDER  TUNNEL-WALL-CONSTRICTION  CORRECTIONS 
IN  TWO-DIMENSIONAL  COMPRESSIBLE  FLOW 
By  E,  B.  Klunker  and  Keith  C.  Harder 

SUMMARY 


Solutions  of  the  first-  and  second-order  Prandtl -Busemann 'iteration 
equations  are  obtained  for  the  flow  past  thin,  sharp-nose,  symmetric, 
two-dimensional  bodies  in  closed  channels.  With  the  use  of  these  solu¬ 
tions  an  expression  is  derived  for  the  tunnel-wall  interference.  The 
tunnel-wall  correction  for  a  parabolic-arc  airfoil  is  calculated  to 
indicate  the  effects  of  compressibility,  ratio  of  the  tunnel  height  to 
the  airfoil  chord,  and  airfoil  thickness  coefficient.  It  appears  that, 
for  cases  where  the  tunnel-wall  corrections  are  significant,  both  the 
second-order  effects  and  the  variation  of  the  correction  along  the  chord 
should  be  considered. 


INTRODUCTION 


The  use  of  wind  tunnels  for  solving  many  aerodynamic  problems  in 
the  hi^- subsonic  speed  range  makes  it  desirable  to  reexamine  the  ques¬ 
tion  of  tunnel- wall  corrections.  The  general  problem  of  tunnel-wall 
interference  in  incon^ressible-flow  fields  has  been  treated  in  refer¬ 
ence  1.  There  the  various  types  of  interference  are  described  and 
equations  are  presented  which  permit  the  correction  of  wind-tunnel  data 
to  free-fli^t  conditions.  These  results  may  be  carried  over  to  include 
the  first-order  effects  of  con^ressibility  with  the  aid  of  the  familiar 
Prandtl-Glauert  rule.  (See  reference  2,  for  example.)  It  is  well-known, 
however,  that  the  linearized  compressible- flow  relations  in  general  do 
not  describe  the  phenomena  accurately  at  hi^-subsonic  Mach  numbers. 

The  exact  analytical  solution  for  the  coir5)ressible  flow  past  an 
arbitrary  closed  body  has  not  been  formulated;  however,  certain  approxi¬ 
mate  methods  prove  useful  for  studying  flow  problems.  The  Raylei^- 
Janzen  method,  in  which  the  initial  step  is  the  complete  incompressible- 
flow  solution,  has  been  used  frequently  for  obtaining  approximate  solu¬ 
tions  to  aerodynamic  problems.  In  reference  3  the  flow  past  a  circular 
cylinder  in  a  tunnel  has  been  obtained  by  this  procedure.  The  Raylei^- 
Janzen  method  is  best  suited,  however,  for  obtaining  flows  past  relatively 
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blunt  profiles  at  low  Mach  numbers;  hence^  this  method  is  not  suitable 
for  evaluating  the  tunnel-wall  interference  for  thin  profiles  at  hi^- 
subsonic  Mach  numbers.  Solutions  of  flow  problems  by  relaxation  methods 
have  proved  useful  in  many  instances;  this  approach  has  been  employed 
in  references  k  and  5  to  evaluate  the  tunnel-wall-constrictlon  correc¬ 
tions  for  the  NACA  0012  airfoil  and  the  Kaplan  bun^,  respectively. 
Although  accurate  solutions  can  be  obtained  by  this  method,  lengthy 
numerical  calculations  are  involved  and  the  result  is  confined  to  the 
specific  profile  and  channel  configuration  treated. 

The  Prandtl-Busemann  small-disturbance  iteration  method  has  been 
found  most  useful  in  the  evaluation  of  the  flow  over  thin  bodies  in  an 
infinite  stream.  This  method  therefore  is  employed  herein  to  evaluate 
the  second-order  tunnel-wall  corrections  for  the  flow  past  thin  sym¬ 
metric  bodies  at  zero  incidence.  The  third-  and  higher-order  solutions 
can  also  be  obtained  by  this  method  although  the  labor  involved  becomes 
prohibitive. 

The  wall- interference  effects  are  often  calculated  from  a  suitable 
system  of  images.  Other  approaches  are  possible,  however,  and  in  many 
instances  prove  more  convenient.  In  particular,  an  integral  representa¬ 
tion  proves  useful  for  evaluating  some  interference  effects.  This  repre¬ 
sentation  is  employed  herein  to  evaluate  the  wall  interference  (con¬ 
striction  effects)  at  subsonic  speeds  for  a  thin,  two-dimensional,  sym¬ 
metric,  sharp-nose  body  in  a  tunnel.  The  method  is  closely  related  to 
the  familiar  source-sink  concepts  for  evaluating  the  flow  over  a  thin 
body  in  an  unbounded  stream.  The  extension  of  the  source-sink  concepts 
to  the  solution  of  the  second-order  Prandtl-Busemann  iteration  equations 
has  been  previously  discussed  in  reference  6.  The  method  of  solution  of 
the  interference  problem  to  the  second  order  thus  parallels  this  work. 

As  an  example  of  the  use  of  the  equations  developed  herein,  the  wall- 
interference  effects  are  evaluated  for  a  symmetric  parabolic-arc  airfoil 
in  a  two- dimens iona],  channel. 


ANALYSIS 


The  exact  nonlinear  equation  governing  the  two-dimensional  flow  of 
a  compressible  fluid  for  the  velocity  potential  is 
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where 

x’,  y’  rectangular  Cartesian  coordinates  in  the  flow  plane 

u%  V*  velocity  conponents  along  the  x*-  and  y’-axes,  respectively 

a  local  speed  of  sound 

and  the  subscripts  x’  and  y*  denote  partial  differentiation  with 
respect  to  these  designated  variables.  With  the  introduction  of  a 
characteristic  length  c/2,  where  c  is  the  chord,  and  the  undisturbed 
stream  velocity  U  as  the  unit  of  velocity,  the  potential  equation  may 
be  written  in  the  nondimens ional  form 


a2/xx  "  ^  V  ^  “  ag/ V  ^ 


,  2x'  2y*  u’ 

where  x  =  — ,  y  =  u  =  v 


20* 

and  0  =  - . 

Uc 


The  Prandtl-Busemann  iteration  equations  are  developed  from  the 
nonlinear  potential  equation  (l)  (see  reference  7^  for  example)  by 
assuming  that  0  may  be  expanded  in  the  form  0  =  x  +  +  02  +  •  •  • 

where  0^+1  Q-nd  its  derivatives  are  small  con5)ared  to  0j^  and  0j^ 

is*  of  the  order  t^  where  t  is  the  thickness  coefficient  of  the  air¬ 
foil.  With  these  assumptions  the  first  two  equations  of  the  Prandtl- 
Busemann  iteration  method  are 


P  ^Ixx  ^lyy  "  ^ 


(2a) 


(2b) 


7  is  the  ratio  of  specific  heats  at  constant  pressure  and  constant 
volume. 
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The  constraint  of  the  tunnel  vails ^  or  the  tunnel- vail  correction^ 
is  expressed  as  the  difference  between  the  velocity  at  any  point  in  the 
flov  in  the  tunnel  and  the  velocity  at  the  same  point  dn  an  unbounded 
stream.  For  aeronautical  applications  this  velocity  increment  is  of 
interest  primarily  on  the  surface  of  the  body,  that  is,  for  correcting 
the  surface  pressures  on  a  body  in  a  wind  tunnel  to  free-stream  condi¬ 
tions.  This  tunnel  correction  on  the  surface  of  the  body  is  found  to 
the  second  order  by  solving  equations  (2)  subject  to  proper  boundary 
conditions, 

Throu^out  this  paper  a  bar  over  a  quantity  (for  example,  ?) 
represents  the  quantity  in  the  channel  whereas  the  same  quantity  with¬ 
out  the  bar  (for  example,  0)  denotes  the  quantity  in  an  unbounded 
stream.  The  differential  equations  (2)  are,  of  course,  the  same  for 
the  flow  in  the  channel  and  in  a  free  stream. 


Boundary  Conditions 


Let  the  equations 


and 


Y  =  tYi(x) 

(-1  <  X  <  1) 

(3a) 

Y  =  0 

(x  ^  1;  X  <  -1) 

(3b) 

define  a  thin  symmetric  body  of  thickness  coefficient  t  lying  on  the 
x-axis  between  -1  and  1.  Then  the  boundary  conditions  for  the  flow 
over  the  body  in  an  unbounded  stresun  are: 

At  infinity 


Ox  ”■  1  }  Oy  —  0  ( 4a ) 

and  on  the  body 

0y(x,Y)  =  y0x(x,Y)  (4b) 

where  Y*  denotes  the  slope  dY/dx.  For  the  flow  over  the  body  defined 
by  equations  (3)  in  a  tunnel  (fig.  l)  whose  walls  are  at  y*  =  ±h’  (or 
y  =  ±h)  the  boundary  conditions  are: 
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At  the  tunnel  wall 


Oy(x,h)  =  0 

(5a) 

on  the  body 

®y(x,Y)  =  Y'Ox(x,Y) 

(5b) 

and  at  x— 

# 

Ox  =  1  >  ^y  =  0 

.  (5c) 

The  boundary  conditions  for  02.  02  order 

respectively,  are: 

t  and  t^ 

At  infinity 

' 

■# 

H 

It 

II 

0 

(6a) 

on  the  body 

01y(x,O)  =  tY* 

(6b) 

and  at  infinity 

. 

II 

n 

0 

(7a) 

on  the  body 

02y(x^O) 

=  tYi'02.x(x,O)  +  32tYi0ixx(x,O) 

y 

(Tb) 

Similarly,  the  boundary  conditions  for  0^  and  02  to  the 
and  t^,  respectively,  are: 

order  t 

As  X— ^±00 

01x  =  01y  =  0 

(8a) 
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at  the  wall  y  =  h 

01y(x,h)  =  0  (8b) 

on  the  body 

01y(x,O)  =  tYi'  (8c) 

and  as 

02x  =  ^2y  =  0  (9a) 

at  the  wall  y  =  h 

02y(^^^)  =  °  (913) 

on  the  body  r 

_  »  ♦ 

02y(x,O)  =  tYi'jZiix(x,0)  +  3^tYi01xx(x,O)  (9c) 

The  solution  of  the  nonhomogeneous  second-order  equation  (2b)  may 
be  expressed  as  the  sum  of  a  particular  integral 

and  a  function  cp2(^^y)  (or  92(^^y))  satisfying  the  homogeneous 
equation 


P^2xx  +  'P2yy  =  ^ 

The  particular  integral  for  equation  (2b)  (reference  8)  is 

i2(x,y)  =  +  |)01  -  I  y01^ 


(10) 


(11) 
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Then  from  equations  (7)  and  (9)  the  boundary  conditions  for  cp2 
92  SLre : 

At  infinity 

'P2X  =  92y  =  0 

on  the  body 

'P2y(^^0)  = 

vhere 

F(x)  s  p2tYi'0ix(x,O)  +  p2tY3^0^^(x^O)  -  +  |)t0i(x,O)0ixy(x,O) 

and  at  x-^doo  - 

92x  =  "^27  =  °  (12a) 

at  the  wall  y  =  h 

9'2y(x,h)  =  f(x)  (12b) 

on  the  body 

■92y(x,0)  =  F(x)  (12c) 

where 

¥(x)  =  32ty^*0ix(x,O)  +  p2tYj^03^^(x,O)  -  M„2^i  +  |)t0i(x,O)  0ixy(x,O) 
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The  tunnel-wall  interference  to  the  second  order  then  is  found  from  the 
solution  of  equations  (2)  subject  to  the  given  boundary  conditions. 


Solution  of  the  Boundary-Value  Problem 

First-order  solution.-  A  solution  of  the  first  of  equations  (2) 
for  satisfying  the  first  two  boundary  conditions  of  equations  (8)^ 

is,  for  y  ^  0, 


0l(x,y)  =  - 


6(0 


cosh  Pa(y  -  h) 
sinh  pah 


cos  a(x  -  I)  d^ 


from  which 


0iy(x,O) 


g(|)cos  a(x  -  I)  d| 


(13) 


The  body  slope  tY^^'  may  be  represented  by  its  Fourier  integral  as 


tYi’(x) 


(l)cos  a(x  -  I)  dl 


(14) 


From  equations  (13)  and  (l4)  the  boundary  condition  for  tangential  flow 
(equation  (8c))  on  the  body  is  satisfied  by  choosing  g(l)  =  ^  Y3^*(l). 
Thus  the  solution  for  Pi{yi,y)  is 


0l(x,y) 


cosh  Pa(y  -  h) 
sinh  pah 


cos  a(x  -  5)  dl 


Inverting  the  order  of  integration  and  integrating  with  respect  to  a 
gives 


0l(x,y)  =  ^  f  Yi'(l)loge  i^jcosh  2X(x  -  I )  -  cos  2XP;0 


y  dl  (15) 
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vhere  X  =  The  solution  0i(x^y)  satisfying  the  boundary  condi- 

23h 

tions  given  by  equations  (6)  may  be  found  from  equation  (15)  by  letting 
h  approach  infinity.  Thus 

0l(x,y)  =  lim  0i(x,y)  =  X  f  Yi*(l)loge  \/(x  -  1)^  +  (3^5) 

h-»oo  P"  ^_1  » 


which  is  the  well-known  solution  for  the  flow  past  a  symmetric  body  in 
an  unbounded  stream  by  a  distribution  of  sources. 

With  Aui  =  0ix(x^O)  -  0ix(x,O),  the  difference  in  velocity  on  the 
body  to  the  first  order  is 

i 

^  J"  Yl’(l)j^  coth  X{x  -  I)  -  — ^dl  (IT) 

Th^  first-order  tunnel-wall  correction  to  the  x-component  of 
velocity  on  the  body  is  given  by  equation  (l?)-  This  solution  was 
obtained  by  considering  the  flow  in  the  upper  half-plane;  the  solution 
in  the  lower  half-plane  is  known  by  symmetry. 

Second- order  solution.-  The  second-order  solution  for  the  flow  over 
a  symmetric  body  by  an  extension  of  the  concepts  of  source-sink  distribu¬ 
tions  has  been  discussed  in  reference  6.  The  second-order  solution  02 
is  given  by  02  =  ^2i  ‘where  cp2  satisfies  the  homogeneous  equa¬ 
tion  (lO)  and  ^2  particular  integral  given  by  equation  (ll). 

The  form  of  the  differential  equation  and  boundary  conditions  for  cp2 
are  the  same  as  for  0j_.  Thus  the  solution  for  02  may  be  written 
immediately  as 


02  =  ^  F(|)loge  ^/(x  -  1)^  +  p^y^  d|  +  +2 

and 

02x(x.O)  =  ^  F(^)  +  0ix^(x,O^ 

(18) 
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(it  should  he  noted  that  the  second-order  solutions  given  here  are  not 
valid  for  hlunt-nose  profiles  since  for  these  cases  the  particular 
integral  introduces  singularities  which  cannot  he  canceled  hy  a  source 
distribution. ) 

The  second-order  solution  for  the  flow  over  a  symmetric  hody  in  a 
channel  may  he  found  in  a  manner  analogous  to  that  for  02*  Consider 
the  expression 


cos  a(x  -  I)  dl  + 


_1^ 

pit 


f(l) 


cosh  pay 


cos  a(x  -  1)  d| 


sinh  Pah 


(19) 


which  is  a  solution  of  the  homogeneous  differential  equation  (lO).  Since, 
from  equation  (19), 


r"  — 

92y(x,0)  =  ij  da J  F(I)cos  a(x  -  I)  dl 


and 


pro  pro 

92y(x,h)  =  diaj  f(|)cos  a(x  -  |) 


d| 


which  equal  F(x)  and  f(x),  respectively,  hy  the  Fourier  integral 
theorem,  the  boundary  conditions  of  equations  (12)  on  the  hody  and  at 
the  wall  are  satisfied.  Inverting  the  order  of  integration  of  equa¬ 
tion  (19)  and  integrating  with  respect  to  a  gives 


cp2 


■— r 

2pit 


F(|)loge 


2p 


from  which  the  boundary  conditions  at 


2\(x  -  I)  -  cos 

2\{x  -  I)  +  cos 
Infinity  may  he 


d|  - 


dl 


shown  to  he  satisfied. 
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The  q_uantity  found  after  some  reduction  as 


02x(x,O)  =  F(|)  coth  X{x  -  I)  d|  - 

r>  00 

f(|)  tanh  Xix  -  I)  d|  +?2x(x,0) 


The  second-order  correction  to  the  velocity  on  the  surface  of 
body  in  a  wind  tunnel  is 

Au  =  Aui  +  Au2  +  0(t3) 

where 


^^2  =  -  02x(^^O) 


or  from  equations  (l?),  (l8)^  and  (20) 

.1 


Au 


coth  \(x  -  I )  - 


1 

X  -  |_ 


dl  + 


F(l)  coth  Xix  -  I)  dl  -  f(|)  tanh  X( 

kl, 


X  -  I)  d5  - 


F(l) 

X  -  I 


*  f)[0l(  x,0)jZ)ixx(x,0)  +  0ix^(x,O)  -  jii(x,O)0ixx(x,O)  - 
^5ix^(x,0)^  +  0(t3) 
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Throu^out  this  paper  the  Cauchy  principal  value  of  all  iii5)roper  integrals 
is  to  he  taken. 


Discussion  of  the  Equations 

The  tunnel-wall- interference  effects  are  often  evaluated  by  con¬ 
sidering  a  suitable  system  of  images.  (See  references  1  and  2,  for 
exan5)le.)  The  first-order  equations  developed  previously  are  consistent 
with  this  approach.  An  expansion  for  coth  X(x  -  |)  (see  reference  9, 
for  exanqple)  is 

coth  h(x  -  I) - .  2  "y  _ » 

k(x  -  I)  n=l  n^it^  +  X,^(x  -  1)^ 

Then,  from  equation  (17) > 


Auj  = 


2kt 

pjr 


-  I) 

k^(x  -  1)^  +  n^jt^ 


d| 


which  is  recognized  as  the  influence  of  source-sink  images  located  at 
y  =  The  series  appearing  in  the  expression  for  Au]^  is  uniformly 

convergent;  hence  the  term-by-term  integration  is  valid. 


The  calculation  of  the  tunnel-wall  interference  for  a  given  profile 
requires  the  evaluation  of  several  Integrals.  The  Integrals  occurring 
in  equation  (21)  are  not  easily  evaluated  for  most  profiles  and  it  is 
preferable  to  expand  the  integrands  in  the  form  of  a  series.  The  series 
given  previously,  however,  is  not  the  most  useful.  The  expansions 


,2n^  2ni 


X  coth  X(x  -  I)  =  .  7"  (.1)-^  .  .^^n-l 

X  -  I  •  (2n): 


tahh  X(.  -  I)  .  (.1)-^  (,  .  I) 

(20).’ 


(A.2(x  -  1)2  <  ^^2) 
2n-li 


y  (22) 


2n-l 


n=l 


(x=(x-l)^<^) 
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(see  reference  10,  for  exaii5)le),  ^ere  the  coefficients  are  the 

Bernoulli  numbers,  are  better  suited  for  calculation  since  the  integra¬ 
tions  are  more  easily  performed.  Moreover  these  alternating  series 
converge  rapidly  and  in  most  cases  only  a  few  terms  need  be  retained. 

With  the  use  of  the  relations  (22)  the  first-order  correction  to 
the  velocity  on  the  body  is  determined  from  equation  (17)  as 


Au- 


=  i  ^  (-1)“-"  f  Yl'dXx  - 1)^-" 

P"  n=l  (2n).'  ^-1 


d| 


(X^ix  -  1)2  <  rt2)  (23) 


or^  after  integration  "by  parts ^ 


n=l  (2n).' 


(A.2(x  -  6)2  <  «2) 

Thus  the  first  term  is  proportional  to  the  area  of  the  body  and  the 
remaining  terms  are  in  a  form  convenient  for  graphical  integration. 

Similarly,  to  the  second  order,  the  correction  to  the  surface 
velocity  is 


00  p2nn  2n-n  pi 

AU  =  ^  21  Yi- 

n=l  (2n).‘  '^-l 


(l)(x  -  |)^-^d|  + 


i  f  F(|')(x  -  |)^"^d|  +  ~  d| 

P"  (2n)!  ^_i  X  -  6 


n=l 


(-1) 


nil  22n(22n  .  i)X^Bn  - 


(2a)} 


f  f(6)(x  -  1)^-^ 
^-1 


d|  + 


+  f)j0l(x,O)0ixx(x^O)  +  01x^(x,O)  -  0i(x,O)lSixx(x,O)  - 

01x2(x,o[|  +  0(t3) 

(x2(x  -  1)2  <y)  (24) 
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The  quantit^s  0^^,  0-]^,  and  their  derivatives  required  to  evaluate  the 

functions  and  f  as  well  as  the  particular  integral  may  he 

found  from  the  equations  (l5)  and  (16),  These  quantities  can  he  expressed 
in  series  form  with  the  use  of  the  relations  (22). 

The  range  of  integration  for  the  integral  containing  f(|)  should 
extend  from  -00  to  ».  The  series  en^ployed,  however,  does  not  converge 
over  this  entire  range.  The  contribution  of  this  integral  is  small 
outside  the  range  -1  to  1;  hence,  with  good  approximation,  the  inte¬ 
gration  may  he  taken  from  -1  to  1.  In  some  cases  the  contribution 
of  the  integral  involving  f(l)  is  so  small  that  it  may  he  neglected 

altogether.  In  this  case,  equation  (2k)  converges  for  “  6  |  <  ’^' 

Since  the  maximum  value  of  x  -  |  is  the  chord,  the  series  form  for  the 
correction  is  applicable  to  cases  where  the  chord  is  less  than  /l  - 
times  the  tunnel  hei^t.  When  the  integral  involving  T(|)  is  considered, 
the  series  solution  is  applicable  to  cases  where  the  chord  is  less  than 

i  \/l  -  Moo^  times  the  tunnel  height. 

The  constriction  effects  for  an  open  tunnel  may  he  found  by  an 
analysis  similar  to  that  given  for  the  closed  tunnel.  Here  the  boundary 
condition  of  zero  normal  velocity  at  the  wall  is  replaced  by  the  con¬ 
dition  of  constant  pressure  along  the  jet  boundary.  To  the  first  order 
the  jet  boundary  is  located  at  y'  =  h'.  However,  for  the  second-order 
solution  the  boundary  condition  must  be  satisfied  on  the  jet  boundary. 

The  location  of  this  boundary  for  the  second-order  solution  may  be 
determined  from  the  first-order  velocities  at  y*  =  h‘ .  The  fact  that 
the  ordinate  of  the  jet  boundary  is  not  a  constant  and  is  initially 
unknown  for  any  iteration  higher  than  the  first  makes  the  solution  of 
this  problem  more  tedious  than  that  for  the  closed  tunnel. 

TUMEL-WALL  CORRECTIONS  FOR  A  SYMMETRIC  PARABOLIC-ARC  AIRFOIL 


As  an  exan^jle  of  the  use  of  the  equations  developed  herein,  the 
tunnel-wall  corrections  are  evaluated  for  a  parabolic-arc  airfoil.  This 
profile  proves  convenient  since  the  first-order  wall  corrections  can  be 
obtained  in  terms  of  tabulated  functions.  Thus  the  rapidity  of  convergence 
of  the  alternating- series  form  of  solution  is  readily  ascertained  by 
comparison  with  this  solution. 

The  parabolic-arc  airfoil  is  defined  by  the  equation 
Y(x)  =  tYi(x)  =  t(l  -  x2) 

=  0 


(-1  <  X  <  1) 
(x  >  1;  X  <  -1) 
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From  equation  (16)  the  u-coii5)onent  of  velocity  on  the  airfoil  in  an 
unbounded  stream  is  given  to  the  first  order  as 


P’f  X  -  5 


2t/o  n  1  -  X 
—  (2  +  X  log  - - 

P«V  1  +  X 


where  (T  denotes  the  Cauchy  principal  value.  From  equation  (15)  the 
■u-coii5)onent  of  velocity  on  the  airfoil  in  the  tunnel  is 


Ut  =  —  Y-,  *  (I )coth  X(x  -  I)  dl  =  -  f  I  coth  X(x  -  l)dl 

Pjf  '^-1  P’*  _1 


and,  after  an  integration  by  parts, 


1  "  P 


(P J  logg  sinh  X(x  -  I)  dl  + 


^l^ogg  |sinh  X(x  -  l)  |  +  logg  jsinh  X.(x  +  l)jj 


Writing  the  hyperbolic  sine  in  terms  of  exponentials,  making  a  change 
of  variable,  and  performing  the  elementary  integrations  gives 


loge  |w| 


a  V  -  1 


dw  +  ^jloge  I  sinh  X(x  -  l)  |  + 


log^  sinh  X(x  +  1)  I  -  2X.X  +  2  log^  ^ 
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where 


w  .  1  - 

a  =  1  - 
b  .  1  - 


The  Spence  integral  is  defined  hy 

X 


El(x)  =  f 

x'  -  1 


Thus  with  the  use  of  this  integral  and  equation  (25)  the  tunnel-wall 
correction  to  the  velocity  is 


-  "‘(‘’O 


2t 

pit 


logg  I  sinh  \(x  -  1)  + 


logg  I sinh  X(x  +  l)|  -  2Xx  +  2  logg  2  -  2  -  x  logg  -  "■  ^ 

1  +  X  I 


(26) 


The  Spence  integral  is  tabulated  for  positive  values  of  x  in 
reference  11.  The  limit  h  is  negative  for  all  values  of  x  for 
X,  ^  0.  However  it  is  easily  shown  that  the  value  of  the  Spence  integral 
for  negative  values  of  the  argument  is  given  hy  the  relation 

RZ(-x)  =  I  RZ(x2)  -  RZ(x)  +  I  RZ(0) 

The  first-order  solution  for  the  tunnel-wall  correction  in  series 
form  is  found  from  equation  (26)  as 


■  -  i  ^  r 

P’'  ^  (2n)!  ^-1 


■d| 


■  ^  h  (-!)“■"  ■  ■  ((X  -  l)""(x  *  2a)  -  .  1)^(,  -  2.^ 

(2a)!n(2n+  I)!- 


(27) 
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Similar  series  expressions  for.  their  derivatives  are 

easily  found.  Then,  vith  the  use  of  these  first-order  solutions,  the 
functions  F(l)^  and  F(lj  maybe  evaluated.  The  integral  Involving 
the  function  f(|)  is  very  small  and  may  be  neglected.  The  tunnel- 
wall  correction  to  the  second  order  is  then  found  from  eq[uation  (24). 


RESULTS  AND  DISCUSSION 


The  first-order  constriction  correction  to  the  x-component  of 
velocity  at  the  surface  of  a  symmetric  parabolic-arc  airfoil  in  a  wind 
tunnel  is  shown  in  figures  2  and  3  for  several  ratios  of  tunnel  height 
to  chord.  It  appears,  at  least  for  the  parabolic-arc  airfoil,  that 
the  error  made  in  neglecting  the  variation  of  interference-velocity 
correction  along  the  chord  is  small  provided  the  chord  is  much  less 
than  the  tunnel  height  and  the  Mach  number  is  small.  For  these  cases, 
the  constant  term  of  equation  (23)  gives  a  good  value  for  the  correction. 
However,  these  are  just  the  cases  for  which  the  correction  is  so  small 
that  it  is  often  neglected.  For  cases  where  the  correction  is  signifi¬ 
cant  (large  Mach  numbers  and/or  small  ratios  of  tunnel  height  to  chord), 
the  error  made  in  neglecting  the  variation  of  velocity  correction  along 
the  chord  may  be  10  to  I5  percent.  Comparison  of  the  correction  (fig.  2) 
given  by  equation  (26)  with  that  given  by  equation  (23)  indicates  good 
convergence  of  the  series  solution.  The  proper  value  of  the  correction 
is  zero  at  a  stagnation  point.  The  fact  that  the  correction  does  not 
approach  zero  at  the  stagnation  points  gives  another  example  of  the 
inability  of  this  approximate  theory  to  represent  adequately  the  flow 
in  these  regions. 

A  series  for  the  first-order  correction  to  the  velocity  is  given 
by  equation  (23);  the  first  term  gives  the  often-quoted  result  that  the 
tunnel-wall  correction  varies  as  1/3^,  whereas  succeeding  terms  contain 
other  powers  of  p.  The  constriction  correction  to  the  second- order  for 
the  parabolic-arc  airfoil  at  midchord  is  shown  in  figure  3*  It  is  of 
interest  to  note  that  the  correction,  to  the  second  order,  at  the  mid¬ 
chord  varies  with  Mach  number  approximately  as  l/P^^  the  second-order 
correction  being  about  20  percent  higher  for  t  =  0.1  and  about  40  per¬ 
cent  hi^er  for  t  =  0.2.  It  therefore  appears  that  the  second-order 
correction  should  be  considered  for  test  conditions  where  the  first- 
order  correction  is  significant. 

The  constriction  corrections  commonly  ejmployed  to  correct  tunnel 
force  or  pressure  data  are  applied  to  the  stream  velocity  and  hence  to 
the  dynamic  pressure,  Mach  number,  and  so  forth  since  the  correction  has 
a  constant  value  over  the  chord.’  On  the  other  hand,  the  corrections 
developed  herein  vary  over  the  chord  of  the  airfoil  and  consequently 
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they  must  he  employed  locally  to  correct  the  surface  velocity  or  pressures. 
These  second-order  constriction  corrections  are  more  accurate  than  those 
normally  employed  hut  their  calculation  and  use  entails  more  labor. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Fields  Va. ^  February  2^  1951 
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Pirst-order  velocity  correction, 


-  Equation  (26) 

- First  term  of  equation  (23) 

-  pipgt  two  terms  of  equation  (23) 


Figure  2.-  Concluded. 
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Figure  3 
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Stream  Mach  number, 


-  Variation  of  first-  and  second-order  velocity  correction  at 
the  midchord  with  Mach  number. 
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Flow.  Subsonic  1.1. 2.1  Wing-Section  Theory 


(Abstract  on  Reverse  Side) 
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